Matrix Element and Strong Electron Correlation Effects in ARPES from Cuprates 
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We discuss selected results from our recent work concerning the ARPES (angle-resolved photoe- 
mission) spectra from the cuprates. Our focus is on developing an understanding of the effects of 
the ARPES matrix element and those of strong electron correlations in analyzing photointensities. 
With simulations on Bi2Sr2CaCu208+i (Bi2212), we show that the ARPES matrix element pos- 
sesses remarkable selectivity properties, such that by tuning the photon energy and polarization, 
emission from the bonding or the antibonding states can be enhanced. Moreover, at low photon 
energies (below 25 eV), the Fermi surface (FS) emission is dominated by transitions from just the 
O-atoms in the Cu02 planes. In connection with strong correlation effects, we consider the evo- 
lution with doping of the FS of Nd2-a:Cea;Cu04±5 (NCCO) in terms of the t-t'-U Hubbard model 
Hamiltonian. We thus delineate how the FS evolves on electron doping from the insulating state in 
NCCO. The Mott pseudogap is found to collapse around optimal doping suggesting the existence 
of an associated quantum critical point. 



I. INTRODUCTION 

Angle resolved photoemission spectroscopy (ARPES) 
has proven to be a powerful technique for probing the 
properties of the superconducting cuprates. In analyzing 
the data, it is important to keep in mind the effects of 
the ARPES matrix element as well as those of the strong 
electron correlation3i*SiMi5iSi2i2iili. With this motiva- 
tion, we present in this article some of our recent work 
directed at understanding these effects in the cuprates. 

In connection with matrix element effects, we discuss 
the ARPES intensity from initial states in the vicinity of 
the M{tt, 0) symmetry point in tetragonal Bi2212. For 
this purpose, one-step photointensity calculations for sev- 
eral different initial state energies for both the bonding 
and the antibonding bands over the photon energy range 
of 10-100 eV are presented. These results give insight 
into how the cross-sections for exciting these bands vary 
with the photon energy on the one hand, and the charac- 
ter and energy of the initial state on the other hand. We 
comment on systematics in the theoretical spectra, and 
provide examples of specific energies or energy ranges, 
which are well-suited for highlighting various interest- 
ing spectral features. We also consider the extent to 
which different sites in the unit cell contribute to the 
photointensity. This question is addressed via computa- 
tions of the dipole matrix element for exciting the bond- 
ing or the antibonding portion of the Fermi surface of 
Bi2212 throughout the Brillouin zone for two different 
polarizations of the incident light. Here the results show 
a remarkable selectivity of the spectra with respect to 
the 0-sites for low photon energies. On the whole, our 
computations indicate that by fine tuning photon energy 
and/or polarization, ARPES can help zoom in on spe- 
cific states and/or sites in the cuprates. We note that 
the methodology used here has been described in detail 
in our earlier publications on the cupratesi*2i2iiiSiSiiiii^. 
The crystal potential used in Bi2212 is the same as that 
used in Ref. [2| and, consistent with ARPES data, it does 
not contain Bi-0 pockets around the Af-point. 



Turning to the electron doped cuprates, we discuss 
the properties of the t-t'-U Hubbard model Hamiltonian, 
keeping in mind the recent doping dependent ARPES 
data on NCCO. The analysis is carried out using the 
framework of the mean-field Hartree-Fock (HE) as well as 
that of the self-consistent renormalization (SCR) theory 
to include the effect of fluctuations. The relevant details 
of methodology are given in Refs.HH anddl It should 
be noted that our comparisons between theory and exper- 
iment on NCCO do not include the effects of the ARPES 
matrix element, but are considered to be adequate for 
our purposes of gaining a handle on the overall topology 
of the FS. An effort to include the ARPES matrix ele- 
ment in the presence of strong correlations is in progress 
and will be taken up elsewhere. The experimental band 
dispersions are used to fit the doping dependence of U, 
which is found to be quite similar for the HF and the 
SCR computations. In this way, we determine the evolu- 
tion of the FS and the Mott pseudogap in NCCO. This 
study also gives insight into the stability of the uniformly 
doped phase (with respect to the formation of competing 
nano-scale orders) upon electron vs. hole doping. 



II. SELECTIVITY PROPERTIES OF THE 
ARPES MATRIX ELEMENT IN BI2212 

Figure ^ presents one-step ARPES photointensities 
for tetragonal Bi2212 in the vicinity of the M{tt, 0)- 
symmetry point at four illustrative initial state energies. 
The k|| -values are chosen for each energy so that the an- 
tibonding (A) or the bonding (B) band is excited as in- 
dicated by open circles in the insert. At energies Ei and 
E2, the A-band is electron- like around F, while at E4, 
it is hole- like around X{t:,tt). Referring to energies, E2 
through E4, the six topmost curves in the figure show 
that the ARPES intensity from either the A or the B 
band is more or less similar for small (compared to the 
bilayer splitting) changes in energy. However, the in- 
tensity is seen to vary dramatically with photon energy 
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FIG. 1: Theoretical one-step ARPES photointensities for 
emission from the antibonding (A) and bonding (B) bands in 
the vicinity of the M{n, 0)-point as a function of the photon 
energy hv in Bi2212. Insert shows the A and B bands and four 
initial state energies E1-E4, at which results for intensities are 
shown (open circles mark the relevant ky -values). Incident 
light is assumed to be polarized along the x-axis, defined as 
the direction of the Cu-O bonds. 



and these variations extend up to the 100 eV upper en- 
ergy limit considered. For example, the A-band is quite 
intense in energy ranges around 20 eV, 47 eV and 88 
eV. The relative intensities of the A and B bands are 
seen to differ greatly. For example, A is far more in- 
tense than B at around 47 eV and 88 eV, while the case 
is opposite around 95 eV. Such special photon energies 
can help focus on the properties of the A or the B band 
in the ARPES spectra. The theoretically predicted en- 
hancement of the A to B features at 47 eV has been 
exploited recently by Chuang et al.^ to resolve bilayer 
splitting throughout the doping range from overdoped to 
underdoped samples, and to adduce that some coherence 
of electronic states between the Cu02 planes persists in 
Bi2212 even in the underdoped regime. 

The preceding discussion makes it clear that the 
ARPES cross-sections for excitation of states near the 
Fermi energy {Ep) depend strongly on photon energy 
and that these photon energy dependencies differ sub- 
stantially between the A and B type states. It is inter- 



esting as well to consider how these cross-sections vary 
with the initial state energy. A sensible energy scale in 
this connection is that provided by the bilayer splitting, 
Abiiaycr- The supcrconducting energy scale, Agupor, is 
~ 30-50% of Abiiaycr- For this purpose, we compare re- 
sults in Figure 1 at the energy Ei (lowest panel) with 
those at E2 — -E4 (upper panels), where Ei is separated 
from E2 — E4 by about Abiiayor- We see, for example, 
that in going from E4 to Ei, the A-band develops peaks 
at around 38 and 77 eV, while the B-band displays new 
peaks at around 18, 63 and 70 eV. On the other hand, 
some energy ranges are less sensitive to changes in the 
initial state energy, e.g. the region from 60-75 eV for the 
A-band and the regions from 35-45 eV and above 75 eV 
for the B-band. These considerations will be important 
in identifying energy regions, which may be particularly 
suitable for implementing recent proposals^^ for develop- 
ing ARPES as a self-energy spectroscopy, where the role 
of kj| and energy dependencies of the ARPES matrix el- 
ement needs to be minimized. 

The ARPES matrix element also exhibits remarkable 
site-selectivity properties in Bi2212 in that, for low pho- 
ton energies (hv < 25 eV), the intensity for emission 
around Ep is dominated by excitations from just the O- 
sites in the Cu02 planes, even though the relevant ini- 
tial states possess an admixture of Cu and O character. 
Our analysis suggests that at higher photon energies, the 
contribution to intensity from Cu sites increases, and at 
around 40 eV, the Cu and O contributions are roughly 
comparable.^ As discussed elsewhere2i^*i^, insight in this 
regard can be obtained by considering the square of the 
dipole matrix element, i.e. 



\Mf = I < */|A-p|*, > I' 



(1) 



where A is the vector potential of the incident photon 
field, p is the momentum operator, and and \1/ / are the 
bulk crystal wavefunctions of the initial and final states 
involved. Fig. |21 presents the intensity based on Eq. 
for emission from the antibonding as well as the bonding 
portions of the FS for two different polarizations of light. 
The contributions to the total intensity arising from the 
0-sites (denoted by Ocu/thin solid line) and the Cu sites 
(thin dashed line) in the Cu02 bilayers are shown. The 
contributions from other sites (not shown) are not always 
negligible as seen by comparing the total intensity curves 
(thick solid lines) with the sum of the thin curves, a point 
to which we return below. Incidentally, O contribution 
in Figs. 121 a) and ^c) is slightly higher than the total 
intensity curve - this reflects interference effects between 
the contributions from different sites when the absolute 
value of the matrix element is taken in Eq. ^ 

In looking at the results of Fig. [3 it is important 
to bear in mind the effect of symmetry on the dipole 
matrix element. In view of the direction of the polar- 
ization vector, the left hand side panels pick up only 
the part My =< '^f\py\'i'i > of the momentum opera- 
tor, p, while the right hand side panels similarly yield 
Mx =< 'i'f\px\'^i >■ Comparing the total intensities in 
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FIG. 2: ARPES intensity defined as the square of tfie dipole 
matrix element, |M|^ (see Eq. 0, for excitation of initial 
states lying along the FS of Bi2212. Thick solid lines give 
total intensity, thin solid lines give partial contribution from 
just the Octi sites in the Cu02 bilayers; the corresponding 
partial contribution from Cu sites is shown dashed. Photon 
energy is 22 eV. The right hand panels refer to light polarized 
along the Cu-O bond in s-direction (denoted by the double 
arrow of the vector potential A), while in the left hand panels 
the polarization vector is rotated by 90°. (a) and (b) consider 
excitation of states on the antibonding FS [thick solid line in 
the insert to (a)]. The ky values along the FS are defined 
by the angle /3 (inserts in a,c), where /3 = 0° corresponds to 
the X(y) — M-direction, and /3 = ±45° corresponds to the 
y(X) — r direction, (c) and (d) similarly refer to the bond- 
ing FS. All intensities are normalized correctly relative to one 
another. 



(a) with (c), or in (b) with (d), we see thus that the 
squared absolute value of My is smaller than that of 
for the antibonding as well as the bonding FS in the 
vicinity of the M(7r, 0)-point. It is straightforward to 
show that the symmetry of the curves in Fig. [3 around 
/3 = is tied to the mirror symmetry of the underlying 
tetragonal lattice through the Y — M line. The mirror 
symmetry through the Y — Y and Y — X lines, on the 
other hand, can be used to deduce that at /3 = ±45°, 
\M^\ = I My I, so that the intensity at /? = ±45° is the 
same for either horizontal or vertical polarization, as is 
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FIG. 3: Dispersion curves for NCCO within the t-t'-U Hub- 
bard model, showing the evolution of the upper and lower 
Hubbard bands (denoted by UHB and LHB) with doping x. 
Results for the mean field computations are given in the left 
column, and for the mode coupling, self-consistent renormal- 
ization (SCR) theory in the right column. Values of x and U 
are marked in each panel. Thickness of bands denotes their 
spectral weight. Horizontal lines mark the Fermi energy. 



seen to be the case with reference to various panels of 
Fig. 12 

Fig. 121 shows that the Cu contribution to the total in- 
tensity is quite small at ^ 22 eV photon energy for either 
polarization for both the bonding and the antibonding FS 
states throughout the Brillouin zone. The Ocu contribu- 
tion is seen to be substantial in all four cases. This is par- 
ticularly the case for polarization along the y-direction in 
panels (a) and (c), suggesting that in the vicinity of the 
M(7r, 0)-point the Ay polarized light would be well-suited 
for probing the character of Ocu electrons. The results 
for the yla;-polarized light, in the right hand side panels 
(b) and (d) show a somewhat more complicated behavior 
in that significant contribution from other atoms (than 
Ocu) in the unit cell is involved. Our analysis indicates 
that these additional contributions arise from O atoms 
in the Bi-0 and Sr-0 layers in the structure. Despite 
some uncertainty with respect to the extent to which Bi 
contributes to emission from the Ep in Bi2212, it seems 
then that polarization will likely be more sensitive for 
investigating the properties of apical 0-sites in the Sr-0 
planes. 



III. MOTT GAP COLLAPSE IN NCCO 

We now discuss briefly the issue of evolution of the elec- 
tronic structure and FS of NCCO with doping and the 
related collapse of the Mott gap around optimal doping 
with reference to Figs. 3 and 4. As already noted, these 
results are based on the t-t'-U Hubbard model Hamil- 
tonian, where the value of U has been used as a fitting 
parameter to reproduce the experimentally observed dop- 
ing dependence of the FS in NCCO observed via ARPES 
experiments. 

The mean field^ and the SCR^ results of Fig. are 
seen to be quite similar, although the specific values of 
U needed in the two cases are somewhat different. The 
SCR computations are of course more satisfactory since 
the effects of fluctuations are incorporated. The mean 
field theory suffers from the well-known problem that it 
predicts the Neel temperature to be too high. However, 
when fluctuations are included (in the mode coupling ap- 
proach), true long-range order only appears at T = 0, 
consistent with the Mermin- Wagner theorem, while a 
pseudogap opens up close to the mean-field Neel tran- 
sition due to short-range antiferromagnetic order. The 
undoped system (x = 0) displays a wide gap between the 
UHB and the LHB. At intermediate doping {x — 0.10), 
the band gap is substantially reduced and electron pock- 
ets appear around the (tt, 0)-point in the UHB. At op- 
timal doping {x = 0.15), the band gap has essentially 
disappeared, and recalling that thin lines possess little 
spectral weight in Fig. 3, the band dispersions resem- 
ble those for the uncorrelated case with a large hole FS 
sheet centered around the (tt, 7r)-point. It is important 
to observe that our analysis indicates that the Hubbard 
U must decrease significantly with doping - by about a 
factor of two from the value of 6< in the undoped case 
to around 3t in the optimally doped systemic. Inter- 
estingly, our finding that the Mott pseudogap collapses 
around optimal doping, where the staggered magnetiza- 
tion goes to zero and the T = Neel ordered state ter- 
minates, suggests the presence of an associated quantum 
critical point (QCP) in the spectrum. The manifestation 
of optimal superconductivity in the vicinity of this QCP 
is consistent with results in other magnetic materialsii, 
including hole-doped cupratesAS. 

Insight into the preceding discussion is provided by 
the illustrative example of Figure 4, which compares 
computed FS maps for the mean field (top right) and 
SCR (bottom left) calculations, with the corresponding 
experimental^^ ARPES intensities (top left) at interme- 
diate doping X = 0.10. In the mean-field calculations, a 
small second-neighbor hopping t" has been included to 
better fit the shape of the FS near the nodal point, and 
the data are broadened to refiect experimental resolution. 
We emphasize that the effects of the ARPES matrix ele- 
ment, discussed in the preceding section, are not included 
in these computations. Nevertheless, we expect accord 
between theory and experiment at the level of providing 
the correct topology of the FS. This indeed is seen to be 
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FIG. 4: Fermi surface of NCCO at a; = 0.10. Experimental 
ARPES results for emission from the Fermi energy in (a) are 
compared with mean field (b), and SCR (c) computations. 
Frame (d) shows a schematic uncorrelated FS (thick solid line) 
and the related "shadow" FS (dashed line) obtained by folding 
the solid line around the magnetic zone boundary. 



the case. Recalling Fig. 3, both theory and experiment 
display the electron pocket around (tt, 0) [with a symmet- 
rically placed pocket located at (0, tt)], arising from the 
doping of the UHB. Additionally, as the Mott pseudo- 
gap decreases, a second hole-like Fermi surface appears 
near {tt/2, 7r/2) coming from the LHB. Notably, evidence 
for two-band conduction has been found in this doping 
range is* 



IV. SUMMARY AND CONCLUSIONS 

In this article, we have briefly discussed two aspects 
of our recent work aimed at understanding the high- Tc 
cuprates via ARPES experiments: (i) Matrix element 
effects in ARPES spectra of Bi2212, and (ii) Evolution 
of electronic structure and FS of NCCO with electron 
doping. 

In considering matrix element effects, we discuss two 
distinct types of selectivity properties of the ARPES ma- 
trix element. Firstly, depending on the character of the 
initial state involved, the ARPES matrix element can 
yield widely differing cross-sections as a function of the 
energy and polarization of the incident photons. This 
effect can be used to discriminate, for example, between 
the bonding and antibonding pieces of the FS in Bi2212. 
Such a theoretically predicted enhancement of antibond- 
ing to bonding intensity at 47 eV has in fact been ex- 
ploited in Ref. 6 to adduce that the bilayer splitting in 
Bi2212 is surprisingly insensitive to doping and that some 
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coherence of electronic states persists across the Cu02 bi- 
layers even in the underdoped regime. Secondly, we show 
that the ARPES matrix element possesses remarkable 
"site-selectivity" in that the emission in the low photon 
energy range (5-25 eV) in Bi2212 is dominated by exci- 
tations from just the 0-sites in Cu02 planes throughout 
the {kx, ky)-j>lane. The polarization vector of the incident 
light may provide some sensitivity to 0-sites in Bi-0 and 
Sr-0 planes. At higher photon energies (above 25 eV), 
Cu sites begin to progressively contribute to the ARPES 
intensity. These results suggest that ARPES could po- 
tentially allow one to focus on the properties of specific 
groups of electrons associated with various sites in the 
unit cell. 

With regard to the electron doped cuprates, we have 
analyzed the evolution of the electronic structure and 
the FS of NCCO with doping within the framework of 
the t-t'-U Hubbard model Hamiltonian in the light of 
related ARPES experiments. Computations have been 
carried out both within the mean field Hartree Fock as 
well as the self-consistent renormalization frameworks. 
We show how the FS develops when electrons are added 
in the undoped insulating state: At first the electrons 
enter the upper Hubbard band around the (0, tt) and 
(tt, 0) points, but on further doping electrons also enter 
the lower Hubbard band first around the (7r/2,7r/2) 



point. Finally, around optimal doping, the FS essentially 
resembles that of the uncorrelated state with a large 
(tt, tt) centered hole sheet. Notably, we find that the 
effective Hubbard U must decrease substantially with 
doping, such that the Mott pseudogap collapses around 
optimal doping, suggesting the existence of an associated 
quantum critical point in the spectrum. Our study also 
indicates that a uniformly doped antiferromagnetic 
state in the cuprates is likely more easily accessible via 
electron doping rather than through hole-doping, since 
in the latter case the system is unstable towards various 
competing orders (nano-scale phase separations). 
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